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ABSTRACT 

In light of the latest IceCube data, we discuss the implications of the cosmic ray (CR) 
energy input from hypernovae and supernovae into the Universe, and their propagation 
in the hosting galaxies and galaxy clusters or groups. The magnetic confinement of 
CRs in these environments may lead to efficient neutrino production via pp collisions, 
resulting in a diffuse neutrino spectrum extending from PeV down to 10 TeV energies, 
with a spectrum and flux level compatible with that recently reported by IceCube. If 
the diffuse 10 TeV neutrino background largely comes from such CR reservoirs, the 
corresponding diffuse 7-ray background should be compatible with the recent Fermi 
data. In this scenario, the CR energy input from hypernovae should be dominant over 
that of supernovae, implying that the starburst scenario does not work if the supernova 
energy budget is a factor of two larger than the hypernova energy budget. Thus, this 
strong case scenario can be supported or ruled out in the near future. 

Subject headings: cosmic rays — galaxies: clusters: general — neutrinos — supernovae: 
general 


1. INTRODUCTION 
The detection of PeV a nd sub-PeV astrophysi- 


cal neutrinos bv IceCube ( 

Aartsen et al. 2013 

Ice- 

Cube Collaboration 

2013 

Aartsen et al. 

m 

4b), 


which more recently has been extended down 


to 10 TeV energies (|Aartsen et al. 2015|) 
major development. 


IS a 

The origin or these neu- 


reviev 

Star-f 

promi 

^s, see e.g., Meszaros 2014 Murase 

2014). 

arming galaxies, especiallv starburs 

ts, are 

sing candidates (Loeb & Waxman 

2006 

Murase et al.||20131 |Katz et al.| 201311'i'amborra 

et al. 

20141 Chang & Wang||2014[ Anchordoqui 

et al. 

2014), in which the major contributing 


novae (HNe) ( Sveshnikova||2003 

Fox et al. 1120131 

He et al. 

120131 liVlurase et al.| 

20131 |Liu et al. 

2014 lE 

lers & Murase 20141) 

or active gaiac- 

tic nuclei 

[ACN) (jMurase et al. 

2014 

Tamborra 


et al. 2014). Other possible candidates include 


low-luminosity classes of 7-ray bursts (Murase 
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2013 Liu fc Wang 2013|), radio-loud ac- 


2013|_witl^ccretion shocks (Keshet et a. 


Muniz & Meszaros 

120041 iVlurase et al. 


Dermer et al. 2014 

Kalashev & 'Troitskv 



ms 


Tum 


Inoue et al. 2005) that may accelerate cosmic 


rays (CRs) to higher energ ies (Hillas||1984 Bier- 
mann & Strittmatter 1987) or oth er OR sources 
such as galaxy mergers m clusters (Kashiyama & 
Meszarc)s||2014 ). 

In this work, we will concentrate mainly on the 
HN/SN origin of neutrinos from pp interactions 
in the starburst and normal star-forming intra- 
galactic material and the intracluster medium. In 
particular, we discuss their implications in light 
of the latest IceCube data in the 10 TeV range. 
Constraints from the diffuse 7-ray background 
measured by Fermi are even more pronounced in 
this case. 

2. HYPERNOVA AND SUPERNOVA 
ENERGY INPUT RATE 

HNe typically belong to a sub-class of broad¬ 
line Type Ib/c SNe with ejecta kinetic ener¬ 
gies of order = 10^^Afc_52 erg, representing 
a fraction ^hn — 4 x 10“^^;m__i.4 of all core- 


Smith et al 


20101 [Guetta fc Della Valle] 


collapse SNe, with substantial uncertainties (|Ar- 
cavi et alj_ 

ik 


2011). The rate of all core-collapse 
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SNe i s 1.06 x 10 ^Mpc ^yr ^ (e.g., 


20141, which implies a local HN rai 
4 X 10-6e^„,_ 


i. 4 Mpc“^yr“^. If the fraction of 
HN remnant kinetic energy transferred to CRs is 
Ecr,hn — 2.8 X 10®^ ergs, the CR energy density 
input rate in the Universe is Ucr — Ecr,hrJi.hn — 
1.2 X ihn-iAEcr,biA erg Mpc“^yr“h Fur¬ 
thermore, if the CRs are protons with a power 
law distribution N{ep) oc £~^ between Ep^min ^ 


Taylor et al. 
e of I^hn ~ 


1 GeV and Sp^max ^ lO^^Eppr eV, the energy 
density input rate per logarithmic interval of en¬ 
ergy is SpQe^ ~ ifcr/C erg Mpc“^yr“\ where the 
bolometric correction is C = ^^{ep^max/sp.min) 
for a spectral index F = ^ Assuming F ^ 2 and 
taking C ^ ISCis the local ( 2 : = 0) CR energy in¬ 
put per logarithmic interval in the Universe due 
to hypernovae is 


— 6.4 X 10 .^cr,/in,51.4 

erg Mpc"^yr"\ (1) 


which is larger than the typical value expected 
for 7 -ray bursts. Conventional SNe will also con¬ 
tribute significantly to lower-energy CRs, hav¬ 
ing a smaller kinetic energy input Ek^sn — 
lO^^Eksn., 51 erg but a larger rate TZsn- The typi¬ 
cal CR energy of SNe is uncertain and could be 
less (e.g. Ecr,sn = 4.8 x 10"^® erg). In general, the 
energy injection rate for SNe is given by 


{SpQe ) = {epQe ) . 

V P^^p/sn C, r PP , ^ P^^pJhn 

<,nn ^sn ^crMn 


( 2 ) 

It is believed that SNe can typically accelerate 
CRs to a maximum energy Sp^max 10^® eV 
resulting in Csn ~ 13.8. The energy input due 
to conventional SNe would be typically larger 
than that of HNe at lower energies. However, 
with the parameters given above, Eq. ® im¬ 
plies the energy input rate of SNe for CRs with 
£p ^ £p,max,sn IS roughly half that of HNe. Be¬ 
low, we leave this ratio as a free parameter. 


3. SHOCK ACCELERATION 

A typical Type Ibc SN has a bulk ejecta 
mass of Mej ~ SMej^o.sAf©, and a HN has 
an average velocity I3ej = {Vejjc) = 6.1 x 

52 ^eya 5 - The postshock random 
magnetic field strength is expected to be ampli¬ 
fied to a fraction es of the postshock thermal en¬ 
ergy, Bs ~ (167reBngTOArc^/3T)^/^, with Ug being 
the interstellar particle number density. The up¬ 
stream magnetic field should also be amplified by 


t It would bed = [1 - (ep,mox/£p,™in)^ *"] £p,mi„(r- 

2)-i if r ^ 2. 


e.g., CR-streaming instabilities, but in any case 
the stronger magneti c fields in starburst galaxies 
may also be enough (Murase et al.||2013). Diffu¬ 
sive shock acceleration in the blast wave leads to 
a power law spectrum distribution N{ep) oc 
typically with F ^ 2 , up to a maximum energy 
~ {i/20)ZeBsRdecPej ( |Drury||1983[ ), or 


-p,max 


£p,hn,max ^ ( 3 ) 

for CRs with charge Z. However, in the follow¬ 
ing work we only consider CR protons. The above 
equation implies that the CRs are accelerated to 
^ 100 PeV as the shock slows and enters the so- 
called Sedov-Taylor phase. The maximum CR 
energy is expected to d ecrease with time during 
this deceleration phase (Drury 20111. Also note 
that while many SNe and HiNe may happen in 
relatively low-density regions such as superbub¬ 
bles, the dependence on Ug is weak. For nor¬ 
mal SNe, using the same parameters except for 
Ek = 10®^ erg, the maximum energy would be 

eV. 


-2/3 


1/6 




-2/3 


£p,sn,max — 1-1 ^ It) 

As we will show below, HNe that occur in star- 
burst galaxies can accelerate the majority of CRs 
which produce detectable high-energy neutrinos. 
Although we use typical numbers for HNe for 
our estimates, in star-forming galaxies hosting 
an AGN the latter may be also contribute as a 


10-100 PeV CR accelerator (|Murase et al.||2014 
Tamborra et al.|[2014 ). 


CRs suffer energy losses both during accel¬ 
eration and after escaping their source. Syn¬ 
chrotron losses are negligible at the energies con¬ 
sidered here, the dominant loss mechanism being 
hadronuclear (pp) collisions. The effective op¬ 
tical depth to pp collisions undergone while ad- 
vected downstream of the blast wave is Tpp^g ~ 
tdyn/tpp ~ KappR{c/V), where k ~ 0.5 is the 
inelasticity and app{£p = 100 PeV) ~ 10“^® cm^ 
(in the numerical calculations presented below, 
we use the energy d ependent inelastic p p cross 
section presented in Kelner et al. 20061. Thus 


1.3 X 10"® A, 


Aff/n r; n„ Vt, which is 


'pp,s x.o ^ /i„,52 0.5 "' 3 . 2.3 

negligible compared to losses during the subse¬ 
quent propagation. Similar considerations apply 
also to the supernovae. 


4. PROPAGATION EFFECTS AND PP 
OPTICAL DEPTH 

The propagation of the CRs in the turbulent 
magnetic field of the host galaxy and galaxy clus¬ 
ter depends, in the diffusion approximation, on 
the strength of the magnetic field B, the CR Lar- 
mor radius r^, and the coherence length of 
the magnetic field fluctuations. At the highest 
energies e^, where /'^(ep) /?> Ic, the CR diffu- 
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sion coefficient is D{ep) cx rL{SpY- At lower 
energies, where r^iep) ^ ic, the diffusion coef¬ 
ficient is D{ep) oc ri(ep)“, where a = 1/3 (1/2) 
for a Komolgor off (Kraichnan) fluctua tion power 
spectrum (e.g., Berezinsky et al.|[l997 |. The two 
regimes can be interpolated as 


UgHappC V[\\iv[td,g,tw,g\ Or 

^ g,2i^g,‘^oBg _^ 'jSp iQ 

'''pp,g ~ 0-55 %,2.3^^g,21^g,20-Bg _3 yEp ]^7 2 

1 ng^2.3Hg^2lVw,3.2 ( 6 ) 


‘pp,g 


BiSp] - (Ep/Ep^:,;) “b (^Sp/Sp^:^^ 


( 4 ) 


where rr iem) = Id^ with U* ~ (l/4)cri(Ep*) 
( Parizot|2004 ). Below we shall use a = 1/3 as an 
example) but the discussion can be generalized to 
a general positive a value. 

After leaving the source (e.g. HNe or SNe), 
the CRs first propagate diffusively through the 
host galaxy or are advected away by a strong 
galactic wind, with typical veloci ties of V,,, 
1500 km/s in starburst galaxies (Strickland & 
Heckman 20091 and V,,, ^ 500 km/s~tor norniaT 
star-t orrning galaxies (Crocker 2012 Keeney et al. 
20061. For a starburst galaxy (SBCj the gas scale 
height Hg ^ 30 — 300 pc may be parameterized as 
Hsbg ~ 300pc ~ 10^^i?2i cm. We assume a mag¬ 
netic field strength of Bg ^ 200 x 10 


in the ranges (Ep > Ep*_g), < £p < £p*,g) and 

(Ep < £w), respectively. 

One can see that starburst galaxies are effi¬ 
cient neutrino factories via pp interactions due 
to their high interstellar gas density. As seen 
in Fig. 1-4, normal star-forming galaxies have 

resulting in 


lower values of £p*^g, £„ 


and 


'pp,g 


-6 


Bg - 3.7 G 


and a coherence length parameterized here as n 


only a modest amount of neutrinos produced at 
high energies. We will show below that if the 
starburst fraction is high and CRs with energies 
up to ~ 10 — 100 PeV are sufficiently confined, 
the majority of the observed high-energy diffuse 
neutrino flux can be explained using HNe or other 
sources in starburst galaxies. 

For the subsequent propagation in the 
galaxy cluster or group, the average mag¬ 
netic field and coherence length are pa¬ 
rameterized as Bri ^ 10~^Bri -fiG and 


^c,g 


10 -Iff, 


30pc ~ 10^°£g,2o cm. Both 
quantities are subject to large uncertainties 
and variations, so that the diffusion coefficient 
adopted here corresponds to the optimistic case. 
For our fiducial starburst galaxy parameters, we 
obtain Ep*^g ~ 1.11 x 10^^£g^2oBg^-3.7 eV and 
D^ g ~ 1.4 X 10^®£g_20 cm^/ s. To ensure GR 
confinement, we require the coherence length t o 
satisfy Ep*_g ^ 10—100 PeV (Murase et al. |2013| ). 
For a normal star-forming galaxy (bb'C), we take 
the typical scale height to be Hg ^ 1000 pc, with 


ru 

Cc,ci ~ 30kpc = lO^^t '23 cm. This im¬ 
plies Ep*,d ~ 5.6 X W^Zi23Bci-3 eV 

and ~ 1-4 x 10^^Z^23 cm^ s”^. For 


a cluster of 10 ^°Mq the virial radius is 
Rci ~ 2.6M35^Mpc ~ 8 X lO^^^M^g^ cm and the 
diffusion time is td^ci = At the maxi¬ 

mum proton energy this is td,ci{£p,m.ax) ~ 2.3 x 


.-1/3 


15 "^23 '^B,-2^k,b2.3^'^ej,l 

s. Similarly to the galactic component 


^p.16.94 

mentioned above. 


there is a spectral break 


_, 

Beck 

2008 

^ ... ........ . 3 .... i 

1 and interstellar medium density of time of CRs ( 

Murase et al. 

2008a 


Assum- 


Given the above, the time for CR diffusive es¬ 
cape from the galaxy can be calculated, which 
for starburst galaxies is td,g = Hg/6Dg ~ 

The time 


Hubble time at the corresponding redshift (i.e. 


t. 


r 2 ^-1/3 


1.5 X i0^^Hg 2ilig,2O^g,-3.7£pd7.2 S- 
for advective escape is t^i^g = Hg/Vw — 6.2 x 
10^^ Hg_ 2 i V// 3,2 s regardless of the CR energy. 
Notice that advective escape dominates diffu¬ 
sive escape from the galaxy for CRs with energy 
Ep ^ £w with 


aged) = /z°° ^ cluster located 

at redshift z = 1 such a break occurs at an 
energy Sp^ci ~ 4 x 10^^^ eV. If the cluster 
break energy is higher than the maximum 
HNe energy, CR diffusion does not significantly 
affect the fraction of CRs that interact in the 
intracluster medium. The cluster pp optical 
depth is again Tpp^ei = ndKappC mm[td,ci,tagej, 


ZeB 


, 1 - 1 / 0 / 


^'11! - 


lOK/gg 
3 c 


1/0 


( 5 ) 


assuming a 
rici ~ IQ-'^ricz -4 cm 
T' — 1) iage ^ id,cldp,‘i^O'X 


typical intracluster gas density 
at high redshifts (e.g. 
so that 


yielding e^, 5.1 x 10^® eV for the fiducial pa¬ 

rameters used here. 

The effective pp optical depth undergone dur¬ 
ing propagation in a starburst galaxy is Tpp g ~ 


7~ppdldp,max) ^ 2.7 X 10 '^rici-dtage/^-S Gyr) 

( 7 ) 

For more nearby clusters Tpp^ci increases as the 
cluster age approaches the local Hubble time, al¬ 
though the density is also redshift-dependent. 
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For a diffusion exponent a different from the 
value 1/3 used as an example above, the values 
of Tpp^g, Tpp^ci,H, etc., are calculated similarly and 
are somewhat different, as can be seen in the nu¬ 
merical results discussed in the next section. 

Further pp collisions occur in the intergalac- 
tic medium after the CRs escape the cluster, 
but with the intergalactic target density riigm = 
2.5 X 10”^(nt,h^/0.022) cm“^, and a total flight 
time limited by tn ~ 10 ^*^yr, the corresponding 
Tpp,igm is negligible compared to the previous two 
contributions. 


5. DIFFUSE NEUTRINO FLUX 

When high-energy CRs undergo pp interactions 
with the ambient intragalactic and intracluster 
medium, charged and neutral pions are created 
which subsequently decay to neutrinos and 7 -rays 
respectively. On average, the resulting neutrino 
and parent CR energies can be related by Si, ^ 
0.03 — 0.05 ep. As a result, the diffuse neutrino 
flux (per flavor per logarithmic interval of energy) 
can be estimated using the CR energy injection 
rate similarly to what is done for GR Bs (Waxman 
&: BahcaII||1997 Murase et al.||201^ , as 






fi,pp 


P 'phys 


6 {l + z’Y 


dt 

dY 


dz\ 


(8) 

where the physical CR energy injection rate per 
energy bandwidth at a given redshift z is related 
to Eq. 0 [ 1 )> cosmological evolution is taken into 
accountoy the scale factor S{z) so the normalized 
physical star formation rate is 

{^pQ£p)phys ^ 


[(epQ, 


^p' hn 


with 


S{z) = 


(1 + z)“^ + 


(SpQe 


1 + z 
B 


(l + zYS{z) (9) 


bp 


l + z 
C 


Ijr] 


( 10 ) 

where a = 3.4, b = —0. 3, c = —3.5, p » 
— lO B 5000, andC 9 (|Hopkins & Beacom 
2006] lYiiksel et al.|[2^ . 


The sum in Eq. Q is over the different galac¬ 
tic and cluster/group contributions. We assume 
a fraction ^sbg of neutrinos are produced in star- 
burst galaxies with the rest ^sfg = 1 ~ Csbg pro¬ 
duced in normal starforming galaxies: 

fpp,sbg —^sbg (l 6 pp.s,«i>s) 

( 11 ) 




Fig. 1. — Diffuse flux per flavor of neutrinos (solid black) 
and 7 -rays (dash-dot) from HNe and SNe, for a diffusion 

coefficient (top): D oc (bottom): D oc in 

both the host galaxy and cluster. For both figures HNe 
and SNe release on average 2.8 x 10^^ and 4.8 x 10^^ ergs 
of CR energy respectively, and the proton spectral index is 
F = 2. The black line with white circ les denotes the mea - 
sured flux of atmospheric neutrinos (|Abbasi et al.|2011||. 
The starburst galaxy scale height, density, and magnetic 
field strength are = SOOpc, risbg = 200 cm~^, and 

Bgbg = 200//G and are represented by red lines. For 
normal star-forming galaxies Hgfg = 1000 pc, Ugfg = 
1 cm~^, and Bgfg = 6 /iG; they are represented by blue 
lines. The contribution from HNe are marked with solid 
lines colored while those from the SNe are dashed. The 
solid green line denotes the total cluster contribution (i.e. 
HNe and SNe from both types of galaxies). Green data 
points correspond to the Fermi measu rements of the ex- 
traga lactic diffuse 7 -ray background ||Ackermann et al.| 
|2014||. Black points correspond to the IceUube measur ^ 
ments of astrophysical neutrinos ||Aartsen et al.||2015|l , 
note that two of the low energy data points are witnm tne 
gray lines of the error bars of the atmospheric flux. 


[e 


sbg I 


+ ^sfg ' 


Note that in the last line of Eq. (11) only CRs 
which escape from the galaxies can contribute to 
the cluster component. 

For our cluster/group parameters and the av¬ 
erage galaxy parameters taken in the diffuse 
neutrino flux per flavor for a D oc diffusion 
coefficient is shown in Fig. [T] (top panel). Here 


/pp.d = (l-e X 
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Fig. 2.— Diffuse flux per flavor of neutrinos (solid black) 
and 7 -rays (dash-dot) from HNe and SNe which release 
on average 7.5 X 10®'^ and 5 X 10^® ergs of CR energy 
respectively, with a phenomenologically motivated diffu¬ 
sion coefficient based on observations of CR diffusion in 
the Milky Way galaxy (see text for details). In this case, 
cluster contributions are dominant at high energies. 


the contributions of the pp interactions in the 
galaxies are indicated both for the supernova and 
hypernova components. In the same figure, the 
resulting diffuse 7 -ray flux is also shown, resulting 
from the corresponding tt® decays and subsequent 
pair cascades in the intergalactic medium, which 
are discussed in ^ A similar calculation for a 
1/2 

D oc Ep Komolgoroff type diffusion coefficient 
is shown in the bottom panel of Fig. 

The situation depends strongly on tlie diffusion 
coefficients of both galaxies and clusters, which 
are uncertain especially at high energies; and, due 
in large part to uncertainties in the magnetic co¬ 
herence length. For example, the diffusion coeffi¬ 
cient for normal galaxies used in Fig. 1 is 10 times 
lower than the value obtained for our Milky Way. 
While this discrepancy is alleviated by inhomo¬ 
geneities, the diffuse Galactic emission suggests 
that the CR spectral break is much lower since 
the observed 7 -ray spectrum is alread y steep at 
GeV energies (Ackermann et al. 2012a|). In Fig. 2, 
we conservatively use the diftusio n coeifficient sug¬ 
gested for our Milky Way by |Blasi fc Amato] 
(2012 1 for normal star-forming galaxies. We_then 
use the scaling relation D oc rL{ep)°‘ (see to 
determine the diffusion coefficient in starbursts. 
Since the break energy is sensitive to the diffu¬ 
sion coefficient, one sees that the diffuse neutrino 
background cannot be explained by star-forming 
galaxies in this case, even with an optimistically 
high fraction of HN kinetic energy converted to 
CRs (i.e. 7.5 x 10^^ erg). At high energies the 
galaxy contribution may not be appreciable. At 
and below ep*_g ~ 1.11 x 10^^ig,2oBg-3.7 eV, 
however, the galactic contribution becomes con¬ 
siderable, Tpp^g(ep*_g) ~ 0.33, overcoming the 




Fig. 3.— Same as Fig. with D oc but top: 

Gba = 0.01 with E^r,hn = 5 x 10®'^ erg, and Ecr,sn = 
2.2x10®® erg bottom: = 0.5 with = 10®^ erg, 

and Ec.r,sn = 2.5 X 10^® erg 

cluster contribution at the same energy. For this 
combination of parameters the cluster and group 
contribution should be dominant, and it is pos¬ 
sible to explain the hard spectrum of the diffuse 
neutrino background. Note that the parameters 
used for the cluster/group contribution to the 
diffuse neutrino flux are optimistic, and massive 
clusters are only a fraction of the total cluster 
population. 

Returning to the parameters used in Fig. 
the flux resulting from average host galaxies with 
a smaller (top) and larger (bottom) fraction of 
GRs produced in starburst galaxies is shown in 
Fig. 1^ Here the fraction of HN/SN GR en¬ 
ergy was adjusted ad-hoc in order to fit the ob¬ 
served neutrino flux with E^r. ,/m = 5 X 10®^ erg 
and Ecr,sn = 2.2 x 10 ®° erg for ^sbg = 0 . 01 , and 
Ecr,hn = 10®^ erg and E^r,sn = 2.5 x 10^^ erg 
for ^sbg = 0.5 respectively. The diffusion coeffi- 

1 /3 

cient was taken to be H oc Ep while leaving the 
remaining parameters unchanged. 

The Fig. and were calculated for “typical” 
star-forming galaxies with parameters as given 
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Fig. 4.— Same as Fig. but for a proton index F = 2.1, 
Ecr,hn = 3.5 X 10®^ erg, and Ecr,sn = 10®° erg. top: 

D oc , bottom: D oc . 

above, and for a proton injection spectrum F = 2 . 
We consider next the SFG and SBG contributions 
using the same parameters, but with a proton in¬ 
jection index F = 2.1, results are shown in Fig. 

H 

The effect of P 7 interactions in the galactic and 
intracluster medium is sub-dominant relative to 
the pp collisions in the relevant energy range, al- 
though it becomes d ominant at very high energies 
( Kotera et al.|[2009 ). 

6 . GAMMA-RAY CASCADES 

The same pp interactions which produce neu¬ 
trinos also produce high-energy 7 -rays with typ¬ 
ical energy e.^ ~ Ep/lO. Note that because of 
this connection, their resulting flux can be re¬ 
lated by When 

7 -rays with energy ^ 100 GeV are injected 
into intergalactic space sufficiently far from Earth 
(i.e. ~ few Mpc), they undergo 77 interactions 
with extragalactic background light (EBL) pho¬ 
tons producing electron/position pairs. The pairs 
scatter additional EBL photons via the inverse 
Compton mechanism generating an electromag¬ 


netic cascade. The resulting cascad ed 7 -ray spec- 


trum takes a universal form, (e.g., 

Berezinskii & 

Sinirnov||l975 

Murase et al. 2U12): 


dN^ 
' de^ 


oc G,., = 


(^) 

(^) 


- 1/2 


l-s 


€7 < £7 


br 


^br 


< £7 < £7“ 


( 12 ) 


and 


given 


^br _ 


The characteristic energies e™' 
above are defined by 1 = [e™‘,z] and 

0.0085(1 -I- ( 0 i^TeV ) respectively and the 

high-energy spectral index is generally taken to 
be s ^ 2. Here is the optical depth for a 
high energy photon traveling through intergalac¬ 
tic spanej_th£ values for which are from model C 
of Finke et al. (2010). 

I'here is also an attenuated component to the 
observed 7 -ray flux from photons with energy 
^ ef,“*which can be calculated similarly to eq. 


_2 ^unatt 
S 7 


47r 


{i+zr 


dz 


o [(l+.s) e., , 2 ] 


2r —2 ^^ \ 

Z_^Ji,pp [SpQsp) 


phys 


(13) 

- e;,,. = 10(l-l-2)£T. 


which combined with Eq. (12) can be compared 
with Fermi-hAT measurement s of the extragalac¬ 
tic d iffuse 7 -ray background (|Ackermann et al. 
'20141. Figures through show our calculated 
diffuse flux of neutrinos and 7 -rays along with 
data from IceGube and Fermi. 


7. DISGUSSION AND SUMMARY 

In this work, we discuss the starburst sce¬ 
nario in light of the new 10 TeV neutrino data. 
Although there are large systematic uncertain¬ 
ties involved in removing the atm ospheric muon 
backgr ound at such low energies (Aartsen et al. 
2014a), it may be challenging to explain the 
diffuse neutrino flux in the whole energy range 
with a single power-law component with F ~ 2. 
Adding the SNe contribution enables us to ex¬ 
plain the low-energy data, but we find that con¬ 
straints from t he diffuse 7 -ray bac kground are 
quite stringent (|Murase et al. 2013|). If the GR 
energy input by SiNe is a factor of two larger than 
that by HNe, the diffuse 7 -ray background is vio¬ 
lated. Additional constraints could be placed on 
the cluster contribution by considering the con¬ 
comitant radio emission, although in the strong 
evolution case these limits are weaker than the 
ones imposed by the accretion shock scenario 
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(Zandanel et al. 2014). From the 7 -ray limits, 
we conclude that it the diffuse neutrino back¬ 
ground in the PeV range originates mainly from 
HNe (and their host galaxies), the HN contribu¬ 
tion should be larger than or at least comparable 
to the SN contribution. However, interestingly, 
in cases where the cluster/group contribution is 
mainly responsible for the diffuse neutrino flux, 
it is still possible for the SN contribution to over¬ 
whelm the HN contribution (e.g. the top panel 
of Fig. [^. 

The strong case scenario, where the ~ 10 TeV 
neutrino data are explained by CR reservoirs, has 
an interesting f eature that can be t ested soon. 


Joshi et al. 2014), some extended sources such 


As proposed by Murase et al. (2013), CR reser¬ 
voirs can give a common explanation for both the 
diffuse neutrino and 7 -ray backgrounds. In gen¬ 
eral, the contributions from starbursts and other 
sources to the neutrino flux above 100 TeV re¬ 
sult in subdominant contributions to the diffuse 
7 -ray background. However, as we show above 
any source which contributes significantly to the 
10 TeV diffuse neutrino flux in the pp scenario 
must also account for almost all of the diffuse 7 - 
ray background. It is commonly believed that the 
diffuse 7 -ray b ackground i s dominated by unre - 


solved blazars (Inoue 2014 Di Mauro et al. 2014), 
implying a cornparatively smaller star burst con- 
tribution. Although there are still significant un¬ 


case scenario can be tested by an improved under¬ 
standing or characterization of the diffuse 7 -ray 
background. 

If, for example, it is proven that blazars are 
responsible for ^ 50% of the observed diffuse 
7 -ray background, the starburst contribution to 
the diffuse neutrino background at low energies 
should be small, especially if the CR energy in¬ 
put from SNe is comparable to or larger than that 
from HNe. Specifically, tighter constraints on the 
the unresolved blazar fraction of the diffuse 7 -ray 
background measured by Fermi and possibly the 
High Altitude Water Cherenkov (HAWC) obser¬ 
vatory (BenZvi 2015), as well as the low energy 


{Si, ^ 100 TeV) spectral shape of the astrophys- 
ical neu trino flux as analyzed by [Aartsen et al. 


(2014aI, will impose limits to the corresponding 


fluxes from galaxy clusters/groups. If there is 
little room for the CR reservoirs, other sources 
need to be responsible for the low-energy neu¬ 
trino component. For instance, there might be 
a significant contribution from Galactic sources. 
Although the Galactic diffuse emission by CRs 
propagating in our Mil ky Way cannot provide 
the main contribution (Ahlers & Murase 2014 


as the I'ernii Hubbles (Ahlers & Murase 2014 
Lunardini et aG 2014 ) or nearby HIN rernnants 


)l'ox et al.|2013'[Ahlers & Murase|2014) could be 
viable. Alternatively, the diffuse neutrino back¬ 
ground might be produced mainl y by hidden neu¬ 
trino sources via jry processes ([Murase fc loka 


2013 


[Kalashev &: Troitsky|[^14[ [Kimura et ar 


2014). The advantage of the strong case consid¬ 
ered here is that it can be tested by multimessen¬ 
ger approaches. It has been commonly believed 
that Galactic GRs come from SN remnants. If 
the diffuse neutrino background is dominated by 
star-forming galaxies, our results imply that even 
Galactic CRs may include significant contribu¬ 
tions from past HN remnants. 

CR acceleration to energies ^ 10^® — 10^^ eV 
has also been proposed in other accelerators 
such as shocks in AGN jets. 


e.f 


Biermann & Strittmatterl 


i2014|), AGN windsl Murase et al.|2014| Tamborra 

et al.||2014 ), and AGJN cores (jStecker et al.||199H 


Hillas 1984 


Uermer et al. 


Ajello et al.J 

2014 

2012 ) and starburst contribu- 

tions (e.g., fAckermann et ai.| 

2012b Tamborra 


et al. 2014), our results imp! 

y that the strong 


Alvarez-Muniz & Meszaros||20U4t |Kimura et al.| 

2014| ). While" n such cases neutrinos can come 
from AGNs themselves, CRs escaping from AGNs 
can also produce neutrinos in intergalactic space, 
which may give a significant contribution to th e 
diffuse neutrino background (Murase et al.|2013). 
Other possibiliti es are accretiW shocks onto clus¬ 
ters of galaxies (IMu rase et al. 2008a Dobardzic 
& Prodanovic 2 galaxy m ergers in clusters 
( [Kashiyama fc Meszaros|2014 ), and 7 -ray bursts 
including tran s-relativistic blNe or low-luminosity 
7 -ray bursts (Murase et al. 2008b Wang et al. 
2007). In principle, the discussion of the clus¬ 
ter/group propagation effects discussed above 
also applies to any intracluster sources. The 
main difference between these other sources and 
SNe/HNe (or sources inside galaxies in general) 
is that the CRs accelerated in the former do not 
undergo pp interactions in the galactic gas, but 
only in the intracluster gas; whereas CRs from 
HNe and galactic sources undergo pp interactions 
in both the host galaxy and the cluster/group. 
This disparity may be relevant at sub-PeV and 
TeV energies, where the spectral shape of the 
neutrino flux can provide clues to the source. In 
this energy range we expect the advective escape 
and Hubble times to dominate the galactic and 
cluster diffusion times respectively, and at differ¬ 
ent critical energies (i.e. Sw and Sp^ci as in (Q. 
Therefore, for galactic sources (that are extra- 
galactic), a soft spectrum is typically expected at 
energies below the maximum acceleration energy, 
with (X (foi' a diffusion time oc e““), 

with a leveling off of the slope to ^constant 
below about ^ 130 (2/(1 -I- z)) TeV (assum¬ 
ing, e.g., a galactic magnetic field strength of 





















































































































200^G and diffusion exponent a = 1/2). Sources 
which release their CRs directly into the intra¬ 
cluster medium on the other hand, are expected 
to produce relatively flat neutrino spectra below a 
break around ^ few (2/(1 -|- z)) PeV, steep¬ 
ening above that. Such a spectral softening of 
the cluster contribution is not seen for the HNe 
model considered here (e.g. in Fig. because 
the break energy due to diffusion is lower than 
the cutoff energy imposed by maximally acceler¬ 
ated CRs. While both and are subject to 
large uncertainties in parameters, the presence of 
their corresponding features could be suggestive 
of sources which are embedded in star-forming or 
starburst galaxies. 


As shown by Murase et al. (20131, in the pp 
scenario, the neutrino spectrum cannot be softer 
than about oc at low energies for the 

corresponding y-ray component not to violate the 
Fermi measurements of the diff use y-ray back¬ 
ground (|Ackermann et al. 2014). At the same 
time, a flat spectrum at moderate to high en¬ 
ergies creates tension with the non-detection of 
neutrinos with energies near the Glashow res¬ 
onance at ^ 6 PeV, which necessitates a 
neutrino spectral sha pe near that energ y steeper 
than £^4)^ oc (Laha et al. 2013). Such a 

spectral break can occur it acceleration stops or 
there is a transition in the diffusive escape time 
(i.e. D{ecr) oc £/), —>• D{ecr) oc £D around 
240 ((1-I-z)/2) PeV (e. g., [Murase et al. 


2013 
in h 


Ar^hordqqui et al. [2014). As can be seen 


also 


ngr~n and [^ the model presented here 
resolve this tension. The neutrino spect 


is flat at low energies £y ^ 130 TeV and soft¬ 
ens to oc £““ slightly before the Glashow 

resonance, while the corresponding diffuse y-ray 
spectrum is below the Fermi measured flux. 

We have shown in (j^that the high-energy dif¬ 
fuse neutrino flux coula potentially be explained 
by HNe, predominantly those located in dense 
starburst galaxies (e.g. the red solid curve in the 
bottom panel of Fig. [h especially for ^ 50 
TeV. For neutrinos wfth this energy and be¬ 
low, the SNe in both starburst and normal star¬ 
forming galaxies contribute significantly to the 
diffuse flux and produce a “bump” in the spec- 
trun^ Reasonable fits by eye are obtained for 
the dirfuse neutrino flux measured by IceCube in¬ 
cluding the latest TeV data by using reasonable 
parameters for the sources as well as the diffu¬ 
sion properties in hosting structures. Such a flux 
also approximates but does not violate the dif¬ 
fuse y-ray background measured by Fermi. This 
does not mean that SNe and HNe are necessar¬ 
ily the only sources contributing to the neutrino 
and y-ray diffuse backgrounds. It supports, how¬ 
ever, the case for GR reservoirs such as clus¬ 
ters and groups being promising sources which 
could contributes at least a significant fraction 
of these backgrounds, wit hout violating both GR 
(Yoshida &: Takarni||2014) and y-ray constraints. 


can 
spectrum 
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